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ABSTRACT 

In the present paper we study generation of the synchrotron emission by means of the feedback 
of Cherenkov drift waves on the particle distribution via the diffusion process. It is shown that 
despite the efficient synchrotron losses the excited Cherenkov drift instability leads to the quasi- 
linear diffusion (QLD), effect of which is balanced by dissipation factors and as a result the pitch 
angles are prevented from damping, maintaining the corresponding synchrotron emission. The 
model is analyzed for a wide range of physical parameters and it is shown that the mechanism of 
QLD guarantees the generation of electromagnetic radiation from soft X-rays up to soft 7-rays, 
strongly correlated with Cherenkov drift emission ranging from IR up to UV energy domains. 

Subject headings: galaxies: active - plasmas - instabilities - radiation mechanisms: non-thermal 



1. Introduction 

Recently an interest to X-ray and 7-ray emis- 
sion from active galactic nuclei (AGNs) has sub- 
stantially increased due to new observational data 
from satellite telescopes. Investigations of cos- 
mic objects passed to a new stage when in 2007 
the AGILE (Astro-rivelatore Gamma a Immagini 
Leggero) satellite and in 2008 the Fermi space- 
craft were launched. These new observational re- 
sults are of fundamental importance in studying 
the emission properties of several X-ray and 7-ray 
sources, like AGNs, pulsars, gamma ray bursts etc. 

Most commonly it is believed that highly rel- 
ativistic electrons do not participate in the syn- 
chrotron emission. Energy losses are so efficient 
that particles very rapidly loose their energies and 
transit to the ground Landau states. In spite of 
this, as i t was shown in a series of articles by 
Osmanovl (l201lll2010allbh:IOsmanov fc Chkhekbe 



(|201l[) ; lOsmanov fc Machabelil (|2010l ). the svn 



chrotron radiation still might account for the ra- 
diation processes produced by very energetic par- 
ticles in AGNs. 



According to iKazbegi et al.1 (|199ll ) in strong 



magnetic fields the plasma may induce the unsta- 
ble cyclo tron waves. On the othe r hand, as it was 
shown bv lLominadze et al. I (|l979h . these waves via 
the QLD affect the particle distribution as along 
as across the magnetic field lines. Such a feed- 
back of cyclotron waves on particles will inevitably 
lead to the creation of the pitch angles, restor- 
ing the synchrotron emission. Apart from AGNs 
the process of the QLD was applied to pulsars as 
well and it was shown that the mentioned mech- 
anism is very efficient for pulsar magnetospheres 
(lLominadze et al.lll979t iMachabeli & Usovl ll979t 



2009 



Malov fe Machabelil 12001 
2010h . ~ 



Machabeli fe Osmanov 



This approach makes possible for the syn- 
chrotron radiation to be a working mechanism 
despite strong energy losses. A very interesting 
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property of the QLD is that it enables to produce 
highly correlated radiation in two different energy 
bands. In this process one has two major forces 
that influence the particle distribution function. 
On the one hand, the diffusion attempts to in- 
crease the values of the pitch angles of resonant 
particles and, and on the other hand, the emitting 
particles are affected by dissipative forces intend- 
ing to decrease their pitch angles. It was shown in 
the papers listed above that under certain condi- 
tions the dissipative and diffusion factors balance 
each other, the physical system reaches station- 
ary state and the pitch angles saturate. Conse- 
quently,the cyclotron emission by means of the 
diffusion redistributes resonant particles and as a 
result, the synchrotron emission is produced. This 
mechanism guarantees strongly correlated radia- 
tion in two different energy bands. 

In the framework of this emission model, the ra- 
diation is generated in two energy bands because 
of the feedback of the excited cyclotron waves on 
particles by means of the QLD, which start to ra- 
diate in the synchrotron regime. In general simi- 
lar emission mechanism can be driven not only via 
the cyclotron waves. Particularly, in the present 
paper, unlike the aforementioned articles, we con- 
sider the feedback of the Cherenkov-drift waves 
on the resonant particles via the QLD. This pro- 
cess should also switch the synchrotron radiation 
mechanism and must provide correlated emission 
in different energy domains. 

The paper is arranged in the following way. In 
section II, we introduce the theory of the quasi- 
linear diffusion. In section III, we apply the model 
to AGNs and in section IV, we summarize our re- 
sults. 

2. Model 

In this section we develop the model of the 
QLD of particles driven by means of the feed- 
back of the Cherenkov drift modes in the mag- 
netospheres of AGNs. Normally, in the men- 
tioned region, with the typical lengthscales I ~ 
10 13_14 cm, values of the Lorentz factors of elec- 
trons may vary from ~ 1 to ~ 10 7 . This range for 
the Lo r entz factors was impl i ed in l Osmanov et aT 



(|2007t) ; iRieeer fc Aharonianl (|2008l) . where where 



pothetical area, where the linear velocity of rota- 
tion exactly equals the speed of light) due to the 
relativistic effects of rotation. We apply methods 
developed in the mentioned papers in the frame- 
work of the present article. For simplicity one can 
assume that the magnetosphere of AGN is com- 
posed of a low energy plasma component with the 
Lorentz factors j p , and a highly relativistic part 
- the beam component with 75. As we have al- 
ready mentioned, the synchrotron radiation is sup- 
pressed for strongly relativistic electrons. In par- 
ticular, one can see that the cooling timescale of 
beam electrons is given by t syn ~ ib'mc 2 / P syn , 
where m is the electron's mass, c is the speed of 
light, P syn s» 2e 4 7 2 _B 2 /3m 2 c 3 is the single par- 
ticle synchrotron emission power, B is the mag- 
netic induction and e is the electron's charge. By 
taking into account typical values of the magnetic 
field in an ambient close to the supermassive black 
hole (SMBH), B w 10 1_4 G, one can see that the 
cooling timescale varies in the following interval 
5 x (10~ 7 — l)s. On the other hand, close to 
the light cylinder zone, due to the efficient cur- 
vature drift insta bility, the magn etic field lines 
significantly twist (|OsmanovM2008| ). therefore, the 
lengthscale is still of the order of I and the cor- 
responding plasma escape timescale, l/c, is of the 
order of 3 x 10 2 ~ 3 s. As we see, the synchrotron 
cooling timescale is much less than the kinematic 
timescale, therefore, in AGN magnetospheres, due 
to the strong energy losses, particles very rapidly 
transit to the ground Landau states resulting in 
the termination of the emission process. 



Ge nerally speaking, as was explained bv lShapakidze et al 
(|2002l ). the necessary condition for the develop- 
ment of the Cherenkov drift instability (ChDI) is 
the presence of the beam component in plasmas. 
Since the magnetic field lines are always curved, 
it is evident that the electrons will drift along a 
direction perpendicular to the plane of the curved 
magnetic field lines with the following velocity 



lb c 



(1) 



it was shown that acceleration of particles becomes 
extremely efficient in the light cylinder zone (a hy- 



where ujb = eB/(mc) is the cyclotron frequency 
of the particle and p is the curvature radius of 
the magnetic field line. It is clear that for highly 
relativistic particles (7?, >> 1) the drift velocity 
becomes significant and the Ch DI with the follow- 
ing resonance condition arises ( Shapakidze et all 
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2002) 



k„v„ 



o. 



(2) 



where fc.. and u are the longitudinal components 
(along the magnetic field line) respectively of the 
wave vector and the particle velocity, and k x is the 
wave vector's component along the drift direction. 

During the process of ChDI the transverse (t) 
waves might be excited , having the growth rate 
(jShapakidze et alj|2002h 



r, 



where 



_ 7TU} blb . 
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7p u x 



(3) 



(4) 



is the frequency of Chcrcnkov emission, ujb = 
\J A^ribe 2 /m is the Langmuir frequency, rib is the 
beam number density, Ak = (k r /k x ) 2 and k r is the 
component of the wave vector along the curvature 
radius of magnetic field lines. Throughout t he pa- 
per we assume Ak = 1 (jKazbegi et alJll9911 ). 

During their motion in a nonuniform mag- 
netic field the charged relativistic particles un- 
dergo two major dissipative forces: H - respon- 
sible for conservation of the adiabatic invariant 



( Landau fc Lifshitd 1 19711 1 



H± = —px, 
P 



H 



PP\\ 



(5) 



and the synchrotron r adiation reaction force 
(iLandau fc Lifshitj[l97l 



P± 



Pa 



a 



2^2 P-L> 



(G) 



where a — 2e 2 u; B /3c 2 , p± and p|| are transversal 
and longitudinal components of the momentum, 
respectively. 

The effect of the aforementioned dissipative 
forces decreases the pitch angles, which inevitably 
causes attenuation of the corresponding syn- 
chrotron emission process. The situation dras- 
tically changes due to diffusion as along as across 
the magnetic field lines, since this process tries 
to increase the pitch angles. Under certain con- 
ditions diffusion and dissipation forces might bal- 
ance each other and as a result the synchrotron 
radiation process will be maintained. The QLD 
influences the particle distribution function and 



the cor responding kinetic equation can be writ- 
ten as dChkheidze et aT1l201lL "~ 
Il979t iMalov fc Machabelife oOl) 



Machabeli fc Usov 
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where / (p) is the distribution function, 
V c / V mc / It 



D 



8tt 



(t) 



2 / e \ 2 W 
\mc/ Tk 



(8) 



are the diffusion coefficients. Here W is the en- 
ergy density of the excited waves and we have 
taken into account th at -Djii = D n = 
(jShapakidze et al.ll2003l ). 

It should be noted that pitch angles ip acquired 
by particles during the QLD are small enough to 
apply the condition, d/dp± » d/dpu, which re- 
duces the kinetic equation (|7|) to the following 
form 



dt 



d 



p± dp A 



(p± 
8 



df° 



p± dpA 



P±D A 



dpi 



(9) 



On the light cylinder lengthscales one can sim- 
plify Eq. © by estimating the ratio F±/H±. Tak- 
ing into account parameters typical for the light 
cylinder area, one obtains 



— - fa 350 x 

F± Aq. 1L40 



10~ 3 rad 
~1> 



10 7 Ric 
x x , (10) 

lb P 



where A 0A = Q/(0.1Sl max ), tt and Vt max ps 
c 3 /(GMbh) are the dimensionless angular veloc- 
ity, the actual angular velocity and the maximum 
angular velocity of the supermassive black hole, 
M BH is its mass, G ps 6.67 x 10 _8 dyne cm 2 /g 2 
is the gravitational constant, L 40 = L/10 40 erg/s 
and L arc the dimensionless and actual values of 
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luminosity of AGN, M 8 = M BH /(10 S M Q ) is its 
dimensionless mass, Mq is the solar mass and 
Ric = c/fl is the light cylinder radius. 

In general, it is evident from Eq. (TTU)) that de- 
pending on physical parameters one has two major 
regimes: (I) F± < H± (II) H x < F±. 

In the framework of the first approximation one 
can neglect its contribution of F± in Eq. (|9]l, 
which for the stationary case (djdt — 0) has the 
following solution 



where 



Cexp 



D±, 



-dp_ 



Ce 



PUT = 



f 2 P D ± ._ 



1/2 



(11) 



(12) 



As we see from Eq. (fTTj) , particles are dis- 
tributed differently for different values of trans- 
verse momentum, therefore, it is natural to exam- 
ine the average value of p± which in turn defines 
the mean value of the pitch angles. A straightfor- 
ward calculation leads to the following expression 



(p 



± I 



(13) 



Jo fi(p±)dp± 

and the corresponding average value of the pitch 
angle 

(P±)l _ 1 Pig 

p\\ p\\ 

For the second case (H± <C F±), the distribution 
function reduces to 



(14) 



fn(P±) = Cexp( ^ 



dp± I = Ce 



where 
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(16) 



and the corresponding mean values of the trans- 
verse momentum and the pitch angle respectively 
are given by 

P±fll(P±)dp± _ v 



(Pa 
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where r(x) is the gamma function. 



(17) 
(18) 



3. Discussion 

As we see from the results of the previous sec- 
tion, the Cherenkov-drift instability strongly in- 
fluences the particle distribution by means of the 
QLD and prevents the synchrotron emission. In 
this section we apply the developed model to 
AGNs and study the production of the nonther- 
mal emission. 

For this purpose we consider the light cylinder 
zone of the magnetosphere. In this area the value 
of the magnetic induction might be estimated in 
the framework of the equipartition approximation. 
Thus, we assume that the magnetic field energy 
density is of the order of plasma energy density, 
which defines the magnetic induction 



B 



lc 



5.5 x Aq.i x LVqG. 



(19) 



As we see the magnetic field is quite strong and 
since in this area particles normally achieve very 
high Lorentz factors, ~ 10 6-7 , the resulting syn- 
chrotron emission should be quite efficient. In par- 
ticular, relativistic electrons moving in a strong 
magnetic field and having the pitch angles ex- 
pressed by E qs. (1141181). will emit the photons 
with energies (|Rvbicki fc Lightman|[l979l ) 
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(21) 



In the framework of the quasi-linear diffusion, 
the problem is usually treated by means of the 
method of iteration. It is evident that the insta- 
bility has to be saturated by means of nonlinear 
effects (the corresponding study is not in the scope 
of the paper) and since there is some energy bud- 
get, it is clear that after the process is relaxed, 
the energy of waves must be of the same order of 
that of plasmas associated wi th the energy bud- 
get (jMalov k, Machabeli|[200lh . Therefore, in the 
framework of the paper we assume W ~ 7f,mc 2 nf,. 

As a first example we consider the case when 
H± exceeds the corresponding component of 
the radiation reaction force. By combining 
Eqs. (|3l8ll2l20p one can study the emission charac- 
teristics. In Fig. [T]we show the dependence of en- 
ergy of radiated synchrotron photons on the beam 
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electron Lorentz factors for different values of 7 P . 
The set of parameters is Ak = 1, Ms — 1, Ao.i = 
1, L = 10 40 erg/s, j p £ [1;3; 10], n b = n GJ , where 
n GJ = QBi c /(2irec) is t he Goldreich- Julian num- 
ber d ensity of electrons (jRuderman fc Sutherland 
19751) . From Eqs. (|12p one obtains that the pitch 
angles vary in the range ~ 3 x 10~ 3 — 2 x 10~ 2 rad, 
that is in a good agreement with the assumption 
H± > Fx (see Eq. HDJ. 

As it is evident from the plots the photon en- 
ergy is a continuously increasing function of 7b, 
which is a natural consequence of the fact that 
more energetic electrons emit more energetic pho- 
tons. ekeV has the similar behaviour with respect 
to the plasma component Lorentz factor. By in- 
creasing 7 P the corresponding photon energy in- 
creases as well. In particular, by combining Eqs. 

(|3I4I8|12[) one can see that ekeV oc D±± ^ lp^- 
As it is clear from the figure, the QLD guarantees 
the emission in the keV energy domain (X-rays), 
which in turn must be strongly correlated with 
the lower energy emission provided by Cherenkov- 
drift instability. From Eqs. (| 1 14|) one can obtain 
energy of photons corresponding to the Cherenkov 
drift waves 



e^«0.06A . l7 - 3 / 2 x 



\ V2 / 



L 



\10 40 erg/s 



1/2 



(22) 

We see that for the mentioned physical param- 
eters the Cherenkov instability produces the IR 
photons with energies ~ (2x 10~ 3 — 6 x 10~ 2 )eV, 
and since the QLD is achieved by means of the 
feedback of these waves on particles, it is evident 
that the IR emission must be strongly connected 
to the soft X-ray radiation (1 — 12)keV. 

Generally speaking, the instability is supposed 
to be efficient if the corresponding growth rate is 
high enough. From Eq. ([3]) we obtain 



T k w 6.5 x 1(T 3 x 7- 1 / 2 x 



lb 



n GJ 



1/2 



( L 



10 7 

1/2 



(23) 



It is clear that the instability timescale, tj nst ~ 
1/Tk is of the order of 500s. On the other hand, 
particles are moving inside the magnetosphere and 
the kinematic timescale (escape timescale) tk in ~ 
Ric/c ~ 5 x 10 3 s exceeds the instability timescale, 




Fig. 1. — Behaviour of energy of radiated syn- 
chrotron photons with the beam electron Lorentz 
factors for different values of plasma component 
Lorentz factors in case of Hx 3> F±. The set 
of parameters is Ak = 1, Ms = 1, Aq.i = 1, 
L = 10 40 erg/s, 7 P £ [1;3; 10], n b = n GJ . 




Fig. 2. — Dependence of energy of synchrotron 
photons on 7b for different values of plasma com- 
ponent Lorentz factors in case of Hx <C F±. The 
set of parameters is Ak — 1, Ms — 1, ^o.i = L 
L = 10 40 erg/s, j p £ [1;3; 10], n b = n GJ . 
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which indicates high efficiency of the Cherenkov- 
drift mechanism. 

From Eq. ([TO)) it is clear that H±/F± is sensi- 
tive with the Lorentz factor of beam components 
and the pitch angles (note that ip itself nontriv- 
ially depends on % see Eqs. (|14I18[) ^ and hence, 
for other parameters the relation between forces 
may be different. Therefore, it is reasonable to 
consider another limit, (Hj_ <C Fx) and see what 
happens in this particular situation. 

In Fig. [2] we show the dependence of energy 
of synchrotron emission on beam Lorentz factors 
for different values of 7 P . The set of parameters 
is A k = 1, M s = 1, A .i = 1, L = 10 40 erg/s, 
7 P G [1;3;10], rib — n GJ . The difference from 
the previous case is jb, which varies from 6 x 10 6 
to 10 7 . From Eq. <jTHJ) it is clear that for the 
mentioned parameters the following approxima- 
tion, Hx -C Fx is satisfied. From the plots we see 
that CkeV lies in the interval ~ (100 — 600)keV, 
which is much higher than the photon energies 
shown in Fig. Q] being a result of a more steep 
dependance CkeV °c lb^ ( see Eqs. (|3l8ll6l21jl ). 
Since the expression of e^y does not depend on 
the beam Lorentz factors, the hard X-ray radia- 
tion is strongly connected with the corresponding 
Cherenkov-drift emission generated in the same 
energy domain as in the previous case: IR band. 

The emission pattern strongly depends on the 
values of the AGN luminosity. In particular, in the 
framework of the model, we use the equipartition 
magnetic field and the induction becomes depen- 
dent on L. Such a dependance is motivated by the 
fact that the diffusion coefficient behaves as, l/£> 8 
(see Eq. ©) which leads to e^ev oc L~ 9 / 8 (see 
Eqs. (1141181201211) ). On the other hand, unlike 
the results shown on the previous two graphs, in 
this case e^y also depend on the AGN luminosity, 
and since synchrotron and Cherenkov-drift emis- 
sion are generated simultaneously, it is interesting 
to investigate them both. 

In Fig. |3]we present the behaviour of the syn- 
chrotron photon energy on the Cherenkov-drift 
emission for different values of angular momen- 
tum. Unlike the previous cases it is worthwhile 
to consider higher luminosity values: [0.5 — 1] x 
10 42 erg/s. We do not examine extreme luminous 
sources, because in this case the Lorentz factors 
of relativistic electrons might be less than consid- 



ered in the present paper. The set of parameters 
is A k = 1, M 8 = 1, A .i e [1;5;10], 7b = 10 7 , 
7 P = 1, rife = n GJ . It is straightforward to show 
that the synchrotron radiation reaction force is 
less than the force responsible for conservation of 
the adiabatic invariant. It is clear from the plots 
that CkeV ( e ev) i s a continuously decreasing func- 
tion, which is a direct result of the fact that with 
bolometric luminosity tkeV decreases (oc L~ 9 ^ 8 ), 
whereas the energies of Cherenkov-drift photons 
behave as L 1 / 2 (see Eq. ([221) ). We also show the 
plots for different values of angular momentum of 
the supermassive black hole. As it is clear, bigger 
the angular momentum, the bigger the resulting 
emission energy for both mechanisms. In particu- 
lar, from Eqs. (|22l8ll6l2ip ) one can see that ek e v 
and e^y both are increasing functions of A. Ac- 
cording to the present results, we see that for the 
mentioned parameters, the quasi-linear diffusion 
provides the generation of synchrotron emission 
from the soft up to hard X-rays: ~ (10 — 250)keV. 
These energies are strongly connected with the 
Cherenkov-drift emission in the energy interval 
~ (1 — 14)eV, corresponding to emission from IR 
up to extreme UV. 

4. Summary 

The main aspects of the present work can be 
summarized as follows: 

1. In the present paper we studied the role of 
Chcrenkov drift emission in maintaining the 
synchrotron regime despite the efficient en- 
ergy losses. 

2. It is shown that in the magnetospheres 
of supermassive black holes the excited 
Chcrenkov drift instability is efficient enough 
to influence the particle distribution by 
means of the quasi-linear diffusion. 

3. In the framework of the model, we examine 
equation governing the process of the QLD. 
It is shown that for physically reasonable pa- 
rameters the synchrotron radiation reaction 
force prevails over other dissipative factors. 
By taking into account this fact the cor- 
responding kinetic equation is analytically 
solved and the emission characteristics are 
estimated 
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4. The emission pattern is investigated for a va- 
riety of physical parameters. We found that 
for typical AGNs the QLD might guarantee 
the excitation of strongly coupled Cherenkov 
drift emission in the eV domain and the syn- 
chrotron radiation in the keV energy band. 

As we have already seen, the Cherenkov-drift 
instability might guarantee the maintenance of the 
synchrotron emission process for a variety of AGN. 
Another important issue we would like to address 
is the problem of radiation spectral index, which, 
in some sense, will complete the problem. It is 
necessary to investigate this task as well, therefore, 
sooner or later we are going to examine it. 
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Fig. 3. — Behaviour of ek e v with e^y for dif- 
ferent values of A in case of H± 3> F±. The 
bolometric luminosity varies in the interval [0.5 — 
1] x 10 42 erg/s. The set of parameters is Ak = 1, 
M 8 = 1, A 0A e [1;5;10], lb = 10 7 , 7p = 1, 
n b = n GJ . 
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